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Abstract. It is shown that unconventional nature of superconducting state of 
PrOs4Sbi2 , a Pr-based heavy electron compound with the filled-Skuttcrudite structure, 
can be explained in a unified way by taking into account the structure of the crystalline- 
electric-field (CEF) level, the shape of the Fermi surface determined by the band 
structure calculation, and a picture of the quasiparticles in f^-configuration with 
magnetically singlet CEF ground state. Possible types of pairing are narrowed down by 
consulting recent experimental results. In particular, the chiral "p"-wave states such 
as Px + ipy is favoured under the magnetic field due to the orbital Zeeman effect, while 
the "p"-wave states with two- fold symmetery such as Px can be stabilized by a feedback 
effect without the magnetic field. It is also discussed that the double superconducting 
transition without the magnetic field is possible due to the spin-orbit coupling of the 
"triplet" Cooper pairs in the chiral state. 
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Recently, the superconductivity has been found in heavy electron compound 
PrOs4Sbi2 with the crystal structure of the filled Skutterudite [Tj. Since the specific heat 
jump AC at the superconducting transition temperature Tc=1.8K is quite enhanced as 
AC/Tc ~ 500 mJ/K^mole, heavy quasiparticles are responsible for the Cooper pair 
formation. Quite recently, a measurement of the longitudinal relaxation rate, 1/Ti, of 
NQR at Sb site has been performed and very unusual temperature (T) dependence was 
revealed both for T < and T > while the normal state properties had been 

known to be also quite unconventional jU Ej. Unconventional behaviors of 1/Ti are 
summarized as follows: 1) Pseudo gap behavior is seen in 1/TiT at < T < 2Tc, in 
which the resistivity p also shows a pseudo-gap behavior PJ. 2) There is no trace of 
the coherence (Hebel-Slichter) peak around T = Tc at all. 3) 1/Ti appears to exhibit 
an exponential T-dependence below 1.3Tc, giving the superconducting gap A in the low 
temperature limit as 2A/A;b^c — 5.3, although a possibility is not ruled out that the 
crossover to the T^-dependence begins to be observed at around T ~ Tc/3 the lowest 
temperature covered by experiments. 

Very recently, the anomaly of specific heat near has been observed, which 
suggests a double transition at T = Td and Tc2 (Tc2 < Td) [SI Ej It also turned out 
very recently on the basis of measurements of the angular dependence of the thermal 
conductivity k under the magnetic field H [5] that there exsist at least two different 
superconducting phases in the T-H phase diagram. In the low-field phase, the 2-fold 
component of along the z-direction is observed as a function of the angle of the 
direction of H around the 2;-axis, while the 4-fold one is observed in the high-field 
phase. The phase boundary approaches the lower critical temperature Tc2 as H ^ 0. 
This is in marked contrast with the case of a heavy electron superconductor CeCoIns 
where the data can be interpreted by a simple "d"-wave model [Hj. These behaviors 
suggest that a novel type of heavy electron superconductivity is realized in PrOs4Sbi2. 

The purpose of this letter is to present a scenario explaining such anomalous 
behaviors in a unified way on the basis of the crystalline-electric-field (CEF) level, 
inferred from the experiment P and theoretical study [Zj , a topology of the Fermi surface 
(FS) offered from the band structure calculations jH], and a picture of the quasiparticles 
of /^-based heavy electrons with a non-Kramers (non-magnetic) doublet CEF ground 
state. 

The CEF level scheme proposed by Bauer et al. [T] in the point group is as 
follows The lowest level is the non-Kramers doublet (P^ in the representation 
of the point group Oh) 



and the first excited state is one of the triplet states P4 or P4 (P4 or P5 for Oh), 




(1) 



r2-3) = 4(l + 2) + l-2)), 



(2) 
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wavefunction of which has the form 

' A?i\-A)-\A)) + 4\\-2)-\2)) 

(3) 

[ B?\ T 3) + 5«| Tl) + ± 1) + 5? I ± 3), 

where the coefficients A^^^'s and B^^^^s are not universal but depend on the details of 
the CEF parameters jZj . The lowest excitation energy of CEF levels has been estimated 
as AcEF — 7K 15 . Other excited CEF levels have excitation energies higher than 
lOOK so that their effects are negligible in the low temperature region of the order of 
Tc. However, the possibility of ri-r4 CEF level scheme, with Acef — 7K, cannot be 
ruled out only from the analysis of the static susceptibility and the specific heat |SI The 
ultrasonic measurements performed quite recently strongly suggest that the is the 
CEF gorund state 

Around T ~10K, these lowest excited CEF levels give considerable contribution 
not only to the thermodynamic quantities, such as the specific heat C and the magnetic 
susceptibility x, but also to the NQR relaxation rate 1/Ti, since the "spin-fiip" process 
can occur among the states forming F^^^ e.g., between | ± 1) and | ±2), giving the NQR 
relaxation. It is noted that each CEF level broadens due to the hybridization with 
the conduction electrons so that the energy conservation law is satisfied in the NQR 
or NMR relaxation process. Indeed, if we assume that the energy level of excited CEF 
level is broadened such that its spectral weight is approximated by the Lorentzian with 
the width AE, and that the processes across the CEF ground states, ((H) and (j2I), and 
r4's are neglected, the imaginary part of the spin-fiip susceptibility lmx±{uj) is given 
simply, in the limit \u!\ T , as 

Imxx(^) ^ const, x -f^^^f^'''''^^^ (4) 

where const, is given by a combinatin of the coefficients A's and -B's in and 
the Clebsch-Gordon coefficeints. Therefore, the NQR/NMR relaxation rate ^ 
Al^Tlmx±{uj) /u due to the excited CEF level is given as 

^ const. X ^-e-^cEP/T (-g^j 



yCEF - — ^^^^ 

The width of CEF level arises from the hybridization between f- and conduction electrons 
and is of the order of the width of the renormalized quasiparticle band. In the present 
case, Acef — 7K is comparable to the bandwidth of heavy electrons, so that AE is also 
expected to be highly renormalized by the correlation effect. 

Namely, if the temperature is decreased well below Acef=7K, the relaxation 
processes are gradually killed, leading to the pseudo-gap behavior (jS)) at such 
temperature region. We have, however, the usual relaxation process due to the 
quasiparticles of the Fermi liquid in addition. Therefore, T-dependence of 1/Ti = 
j^^jiCEF _|_ i^j'^p at T < Tc will be rather complicated one, since both contributions, 
from r4 CEF level and the quasiparticles (1/T[^p), to 1/Ti are decreasing 
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in such T-region with different T-dependence in general. In particular, one has to be 
careful in drawing the structure of the superconducting gap from the T-dependence of 
1/Ti at Tc/4 < T < Tc. 

In this Letter, we discuss the nature of the gap structure specific to the present 
system. First of all, it may be reasonable to assume that the strong on-site repulsion, 
the possible origin of the heavy electron state, cannot be avoided in a manifold of the 
conventional s-wave pairing state. This is also consistent with the absence of the co- 
herence peak in 1/Ti [2\ although the possibility of the strong coupling effect is not 
completely ruled out. Another crucial aspect of PrOs4Sbi2 uncovered by the band struc- 
ture calculation is that the FS of the heavy electron band is missing in the directions 
of [1,0,0] and [1,1,1] and their equivalents as shown in Fig. ^ |HJ. There exists a small 
FS surrounding the F-point whose mass is not heavy [S] and has been detected by the 
de Haas- van Alphen experiment jTTj. Due to this porous structure of FS, even the 
anisotropic pairing state can have finite gap over the FS. However, even with such FS 
features being taken into account, the anisotropy of the gap due to such features of FS 
does not seem to fully explain the exponential-like T-dependence of 1/Ti. 

Odd-parity pairing due to quadrupolar fluctuations 

We adopt here the odd-parity pairing to explain the unconventional nature of 
the superconducting state mentioned above. There are at least three circumstantial 
evidences for favouring the odd-parity pairing. First, the pairing interaction should be 
mainly mediated by the mode which gives rise to mass enhancement of quasiparticles, 
the quadrupolar fluctuations in the present case. Quadrupolar susceptibility XqIq) is 
expected to be enhanced at large wave vector because the main FS has the nesting 
tendency as shown in Fig. |21 which can induce the attraction in both the "d"- and 
"p"-wave channel since the spin factor [a ■ a') does not exist in contrast to the case of 
spin-fluctuation mechanism f2lE]- Second, a scenario for the bouble transition is more 
easily constructed in the odd-parity pairing with degeneracy due to the time-reversal 
symmetry than the even-parity pairings. Third, the so-called Maki parameter K2 under 
the magnetic field exhibits no paramagnetic limitation |14j . 

As mentioned above, the pairing is also expected to be induced by exchanging the 
quadrupolar fluctuations of the non-Kramers doublet F^j. The quadrupolar coupling 
between essentially localized 4f^-states and the quasiparticles containing considerable 
weight of the conduction electrons, arises through the hybridisation with the local 
symmetry of Fg^"* and Fg^"* in the cubic representaion of CEF state of j = 5/2-manifold, 
as discussed by Cox in the context of the quadrupolar Kondo effect [15 . The fact 
that the heavy quasiparticles contains the considerable weight of conduction electrons 
is a salient feature of f^-based heavy electron state, which is in marked contrast with 
the f^- or f^-based ones where the quasiparticles are dominated by f-electrons. It is 
also consistent with the result of band structure calculation which shows that the f- 
component of the heaviest band at FS is only several % |Hj. The propagator of the 
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quadrupolar fluctuations Xq{1j i"^™) may be given as 

XQ{q,iuj„,) = . (6) 

1 - gmo{q, luj^) 

where Xq"^ and Ho denote the propagator for local fluctuations of quadrupolar moment, 
and the polarisation function of quasiparticles, respectively, and g is the coupling 
constant among them. 

As shown in Fig. |2l the FS has a nesting tendency and Ho is expected to have peaks 
at g = (7r/2, 7r/2, 0), and q = (7r/2, 37r/2, 0), and their equivalent positions. It is noted 
that the FS is rather flat in the 2;-direction near the nesting position as can be seen in 
Fig. H Then, the pairing interaction in the static approximation, T{q) ~ 5'^Xq(Q')0)5 
can be parameterised as 

T{q) = Fo - ri[cos(2g^) + cos(2gy)] + T2 cos y cos y 

+ (cyclic permutations of g^, qy, and g^), (7) 

where Fj's are positive constants and r2 is rather smaller than than Fi. The term of F2 
represents the effect that tendency of the nesting at q= (7r/2, 37r/2, 0) is less than that 
at g* = (7r/2, 7r/2, 0) as seen in Fig. |21 By putting q = k ~ k', F(g) is represented near 
the peak as follows: 

T{k-k') = To-Ti[cos2{k,-k'^) + cos2{ky-k'y)] + --- (8) 
This gives the attractive interactions in the following channels: 

"d"-wave; cos(2A;^) — cos(2A;j^), etc., (9) 

"p"-wave; sm{2kx), sin(2A;j,), sin2A;2, (10) 

Among these states, sin(2/ca.) and its equivalents, sin(2A;j^) and sin(2A;2), will be the most 
favourable ones because they have maximum amplitude on the FS. Indeed, other states 
have more nodes on the FS as seen in Fig. |21 

The simplest odd-parity states with "equal-spin-pairing" (ESP) allowed in cubic 
symmetry are given as follows: 

Ak = A[px{k) + epyik) + e^p,ik)]i{ayax), (11) 
Ak = A[p,(fc) + ipy{k)]i{aya,), (12) 

Here, aj is the j-th component of the Pauli matrix, e = e'^'^/'^, and p(/c)'s are bases of 
irreducible representations with "p"-symmetry: Px{k) = a/2 sin (2 fc^;), etc.. These gaps 
vanish along the direction of [1,1,1] or [1,0,0], and its equivalent direction. However, 
since there exists no FS in those directions, 1/Ti exhibits an exponential T-dependence 
in the lowest temperature region in spite of the anisotropic gap. It is noted that the 
Fermi surface of light electrons, detected by the de Haas- van Alphen experiment JT], is 
closed surrounding the F-point and these gaps have nodes at points on this FS. However, 
since 1/Ti is proportional to the square of the density of states at the Fermi level, the 
effect of such light electrons should hardly be seen by the T-dependence of 1 /Ti . 
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Other possible states in the odd-parity manifold are 

Afc = Ap^(A;)i((Ty(Ja;), and its equivalent ones. (13) 

Such states are less favourable compared to the chiral states (|TT|) and (fT^ in the so- 
called weak-coupling case where the feedback effect is not taken into account. This can 
be seen from the structure of the GL free energy. For instance, in the case of odd-parity 
class of ESP with "p"-symmetry, it is given as follows IT7]: 

Fgl(A,, A„ a,) = Fo + $(1 - V<I>){\A,\^ + |A,|2 + |A,n 
+ iXdiag(|A,|^+|A,|^+|A,|^) 
+ Xoff{|A,nAj2 + |AJ2|A,|2 + |A,nA,|2 

+ 2[Re(A,.A;)]2 + 2[Re(A,.A;)]2 + 2[Re(A.A;)]2|, (14) 

where A's are the coefficient of each (normalized) irreducible representations of the gap, 
V is the strength of pairing interaction of "p" -symmetry, and 

» ^ EMk)?'-^^^%^< (16) 

k 

k \ J 

Since Xdiag > Xos due to the Schwarz inequality, the gap (fT^ cannot minimize Fql in 
general. In case Xofr < Ixdiag, the gap (fTTj) minimizes -Fgl, while the gap fT^ minimizes 
Fgl in case Xdiag > Xos > |xdiag- 

The low-field phase, in the T-H phase diagram having 2-fold symmetry is con- 
sistent with the gap ()13|1 . which cannot be realized in the weak-coupling theory. This 
is also the case in the singlet manifold. The so-called BW-like state is known to be the 
most stable state in the weak-coupling approximation, and there is no reason in princi- 
ple to rule out its possibility from the first. However, BW-like state looks inconsistent 
with the thermal conductivity measurement under the magnetic field |3] and other ther- 
modynamic measurements. 

Feedback effect 

In order that the gap with 2-fold symmetry such as !^l'S\i to be realised, we need 
a feedback effect. Among them, the following mechanism may be promising. The 
polarisation function IIq appearing (jH)) in the superconducting state is given as 

U (r, n\ ^ ^ ^k^k+q ~ ^k^k+q + ^k^k+q ( . o\ 

no(g,o) = -l, EE^^^E.E , ^ (18) 

j: k k+q^ k k+qy 



where = y^^S+|AgP. If the nesting were perfect at g = (7r/2, 7r/2, 0), and 
q = (7r/2, 37r/2, 0), and their equivalent positions, the following relations would hold. 
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^fe+g- = -^fc, Afe+g' = -Afc, and E,:^^ = E^:, for k near the FS. Then, the expression, 
flTS|) . would be reduced to 

no(g,o) = ^E ..g^fLw. (19) 

Then, the polarisation mediating the pairing interaction depends on the type of pairing 
itself. Indeed, the pairing (fT!^ is expected to give larger (jTH]) than the pairing (fT^ . 
because the gap function of (fT^ 

|A^|2 oc 2sin2(2A;^) (20) 

vanishes on the planes, = 0, ±vr/2, ±7r, and ±37r/2, which pass through near the FS, 
while the gap function (fT^ 

I A^|2 oc [sin2(2A;^) + sm^{2ky)] (21) 

vanishes only on the lines, {kx,ky) = (±7r/2, ±7r/2), (±7r/2,±7r), (±7r/2, ±37r/2), etc., 
which are located away from the FS. Although the explicit band structure calculation is 
hard in practice for the moment, the tendency mentioned above is expected to remain 
valid. Therefore, the state (fT^ may be stabilized against (fT2|) by the feedback effect. 
In the BW-like state, no(g, 0), (fT^. is suppressed more severely than in the state (fT^ . 
and is destabilized against (fT^. 

This kind of feedback effect is an analogue of that due to the ferromagnetic spin- 
fluctuation mechanism discussed in superfiuid ^He [IHllini, in which the spin-fluctuation 
spectrum depends on the gap structure of the triplet states. 

Double transition due to spin- orbit coupling 

The spin-orbit interaction Hso due to the mutual Coulomb interaction between 
electrons and relative motion is given by 

= E E + ^j) ■ [^L- X [(2^ - m - , (22) 

where /ib is the Bohr magneton, m electron mass, mband the band mass, and g is 
defined as g = /^eff/z^B, A^eff being the effective magnetic moment fies = (6/7)| (j^) I/^b- 
The appearance of the factor m/mband in ()22|1 can follows from the Ward-Pitaevskii- 
identity [20] • By the procedure similar to that described in Ref. |16j for the dipole 
interaction, the interaction (j^^ leads to the spin-orbit free energy Fg^ for Cooper pairs 
which is spin triplet and chiral, such as ()12|1 . with the pair angular momentum M as 
followsj^n]: 

F,, = -gU^dxd*)-i, (23) 



where 



m 20,, _ m 
9d X —{ig - 1) = X { 

"^band -J "^band 



20 37 
— X — , for r 

3 7 ' 



(2). 



(24) 



20 5 
— X -, for r 
3 7' 



(1) 
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where gu is the strength of the dipole couphng in the "ESP" -superconducting state, and 
we have used {jz) = ±11/6 for quasiparticles consisting of T^'^ f^-CEF state, and ±1/2 
for rg'\ The free energy due to the dipole-dipole interaction is given as [E] 

Fj, = -^g^\{d (25) 

Therefore the spin-orbit interaction, in the non- unitary state with |(ix(i*| ~ 1, dominates 
the dipole-dipole interaction, in the unitary state with \d-i\ = 1, because gso far exceeds 
g-D in Fg^^-band considering and rn/m\,s,ad ~ 0(10"^). Following the calculation in 
the spherical model jTHI, gu is given by 

= Y^TJTc ^-^-"'ff^^ [ln(1.14ec/A;BTe)]' , (26) 

where Np is the density of states (DOS) of the quasiparticles, and Fcond is the 
condensation free energy 

where k is the average of square of the magnitude of normalized gap function over the 
FS. The second factor in is estimated as 

3.1filsNF [ln(1.14ec/A;BT,)]' ~ (^] x 1.4 x IQ-^, (28) 

where we have assumed that the renormalized Fermi energy is ep ^ 10^/300 K, the 
number density of quasiparticles N/V = 2/(2/v^rpi._pr)^, Tpr-Pr being the distance 
between two nearest Pr ions. Therefore, the spin-orbit coupling ^fso, is estimated 
as 



^ cond 

'/ ^ 

- T/Tc rriband 



r 20 


37 


X 


— X 1.4 X 10 


3 


7 


20 




X 


- X 1.4 X 10" 


I 3 
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4.9 X 10-\ for F?^; 



6.6 X 10-2 for F^^^ 



The free-energy difference between (fT!?|l and (fT^ . its non-unitary version with |(i x 7^ 
0, is of the order of 10% of -Fcond in general, and m/mband ~ 0(10^^) according to the 
band structure calculation [Hj. Therefore, if the stable state is (fT!?|) due to the feedback 
effect as discused above, there occurs a double transition with splitting of the transition 
temperature being (T^ — Tc2)/Tci — several %, because the state (fT^ of non-unitary 
version is stabilized, due to the spin-orbit interaction (j2Sl), against (fT!?|l which is real 
state and has no spin-orbit coupling such as (j2Sl)- This width of splitting of double 
transition is consistent with the experimental observations [Hllll- Especially, the results 
by Aoki a/ jS], suggesting the double transition remains rather robustly under the 
magnetic field, can be explained by the present mechanism. 

The chiral state (fTT|) has also intrinsic magnetic moment along (1,1,1) direction, so 
that it can also give rise to the double transition as above. However, this state gives 
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the angualr dependence of the thermal conductivity opposite to the observation in the 
high-field phase fS| although the 4-fold behaviour is expected. 

Two phases in T-H phase diagram 

Finally, the multiphase diagram in T-H plane determined by the thermal 
conductivity measurements under the magnetic field j3] may be understood as follows: 
A crucial point is that the triplet state ()12|) has the intrinsic magnetic moment Mjn 
associated with the intrinsic angular momentum Lin as Mm = /UB("^/''7^band)-^in/^, where 
m* is the effective mass of heavy quasiparticles, and Lin is given as 



where A^in is the order of the superfluid electron density j^H 122 12^]; while the 
lively disputes were performed concerning the size of A^in, whether A^in ~ C'(A's) or 
0{Ns ■ (Tc/ep)") with n = 1 or 2, about a quarter of century ago in the context of 
superfiuid ^He Very recently, the reality of this intrinsic magnetic moment has 

caused a renewed interest in the magnetic property of the chiral superconducting state 
of Sr2Ru04 [25 . At low enough temperature T -C T^, the intrinsic magnetic moment is 
Min — (A^/2)/iB(m/mband)^- Therefore, the state (|T^ is stabilized under the magnetic 
field H over the state (fT^ . which has no intrinsic magnetic moment. The transition 
between the two states occurs when 



where the left-hand side represents the difference of the condensation energy between 
the two superconducting phases at T <^ Tc, and the right-hand side the energy gain 
in the chiral state (jl2p with the intrinsic angular mementum in the magnetic field H. 
We have included in the so-called the demagnetisation factor which depends 
on the sample shape. In ()31|). the energy due to the magnetic field arising from the 
intrinsic magnetisation itself is neglected because it is much smaller than the external 
field ~ if in question. The magnetic field giving the phase boundary beween low- 
and high-field phases, determined by the thermal conductivity, roughly agrees with 
the present estimation is in the same order as given by (jSH), because Ap ~ N/ep, 
ksTc/ep ~ 10 X {m/m*), m*/mband ~ 10, and A'd ~ 10^^. A crucial prediction of 
the present scenario is that the phase boundary in the T-H plane is dependent on the 
sample shape through the demagnetisation factor Aj. 

The angular dependence of the thermal conductivity Kz reported in Ref. ^ may 
also be explained by the present scenario. If the state (fT^ . oc Px{k), is reahsed in 
the low-field phase due to the boundary effect, which works to align the extention of 
pair wavefunction, takes maximum (minimum), when the magnetic field B is B \\ x 
{B II y), in consistent with Ref. This can be understood by applying the argumet 
similar to Ref. jSj. If a type of the state ()12|1 is realized in the high-field phase, the 
free energy takes minimum when the quantisation axis of intrinsic angular momentum 
is parallel to B for which Kz is minimum Therefore, when the direction of B is 



'm 




(30) 




(31) 
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Figure 1. Fermi surface of PrOs4Sbi2 relevant to tlie heavy electrons given by band 
structure calculation FS is missing in the direction of [1,1,1] or [1,0,0], and its 
equivalent ones. 

rotated by angle from the x-axis in the plane perpendicular to 2;-axis, Kz increases 
and reaches maximum at = 45° above which the stable quantisation axis changes 
from X- to y-axis, and then Kz decreases up to = 90° where the stable configuration is 
reached again. Namely, the 2- fold dependence of Kz{4>) taking maximum at = 0° can 
be possible in the low-field phase, and the 4-fold one taking minimum at = 0° and 
= 90° in the high-field phase [Sj. 

In conclusion, it is remarked that the magnetic susceptibility, both longitudinal and 
transverse, can be enhanced by electron correlations even if the mass enhancement arises 
from the degeneracy due to the non-Kramers doublet, i.e., electric quadrupolar moment, 
provided that there exists a perturbation which breaks the particle-hole symmetry, such 
as the repulsion among conduction electrons, as shown by the numerical renormalization 
group calculations for the impurity model [26J. This is in marked contrast with the case 
of heavy electrons based on f^-configuration with the singlet CEF ground state |27j . 
where the static susceptibility along the easy axis due to quasiparticles is not enhanced 
while the NMR/NQR relaxation rates given by the dynamical transverse susceptibility 
is enhanced in proportion to a square of the mass-enhancement factor as observed in 

UPta I2H1- 

We have benefited from informative conversations with Y. Aoki, T. Goto, K. Izawa, 
Y. Kitaoka, H. Kotegawa, Y. Matsuda, H. Sato, H. Sugawara. This work was supported 
by the Grant-in-Aid for COE Research Program (No. 10CE2004) by the Ministry of 
Education, Culture, Sports, Science and Technology. 
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Figure 2. Fermi surface sliown in Fig. cut by the k^-ky plane k^ = 0. Tlie nesting 
tendency in the heavy electron band at (k^, ky) = (7r/2, 7r/2) and (37r/2, 7r/2) as shown 
by arrows, and their equivalent positions, remains in the direction rather rubustly for 
finite k^ ^Q. 
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